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characterizing the degree of "hardness" of the instanton product is presented. 
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I. INTRODUCTION 

In the Standard Model, both the strong and 
the electro-weak interactions are described by non- 
Abelian gauge theory. In these theories, the ground 
state has a rich topological structure, associated 
with non-perturbative fluctuations of the gauge field, 
called instantons 0, 0, LH , which represent tunnel- 
ing transition between topologically non-equivalent 
vacua. In the strong sector, described by QCD, 
instantons are non-perturbative fluctuations of the 
gluon field. Deep-inelastic scattering (DIS) offers 
a unique opportunity to discover the processes in- 
duced by QCD instantons. The rate is calculable 
within "instanton-perturbative theory" and is found 
to be sizable |g, |y, [jj • 

QCD-instanton induced process leads to a char- 
acteristic final state, which allows instanton induced 
events to be distinguished from the normal DIS pro- 
cesses. For a leading graph in QCD-instanton in- 
duced e-p collision, cf. Fig. 1, the incident lepton 
emits a photon, with 4-momentum q, which in turn 
transforms into a quark-antiquark pair. One of these 
quarks with 4-momentum q" hadronizes to form the 
current jet. The other quark, with 4-momentum q' , 
fuses with a gluon (4-momentum g) from the proton 
in the presence of an instanton. The q'g interac- 
tion is called the hard subprocess of photon-gluon- 
fusion process. For simplicity, we will in the follow- 
ing call the hadron system produced from the fusion 
of quark-gluon in the presence of an instanton as 
instanton final state (IFS). 

The phenomenological characteristics of instanton 
induced events can be summarized as follows Q: 

• In the hard subprocess exactly one qq pair of 
each of the n r kinematically accessible quark 




FIG. 1: The leading QCD-instanton induced process 
in the DIS regime of e-p scattering. 

gluons plus 2n/ — 1 outgoing (anti-)quarks give 
rise to a high multiplicity final state. 

• The particles produced from the quark-gluon 
fusion in the presence of instanton, i.e. the 
IFS, are expected to be isotropically dis- 
tributed in their center of mass frame. 

• From theoretical point of view, the most 
prominent characteristic of instanton induced 
events, which is also the most difficult feature 
to be observed experimentally, is the chiral- 
ity violation. In instanton induced events, all 
quarks produced in the hard subprocess are 
emitted with the same handedness, that means 

7* + <7 ffi - ^(qR + q^+ngg, (I^I,R^L). 

n t 

(1) 



flavors participates in the hard subprocess of Although the existence of instanton is required by 
each event, either as incoming or as outgoing the Standard Model, the experimental evidence is 
object. This feature is often termed as "flavour still lacking. HERA offers a unique chance to dis- 
democracy" . The number of gluons emitted cover instanton, which would be a confirmation of an 
per event follows a Poisson distribution with essentially non-perturbative Standard Model predic- 
average value (n g ) ~ 0(1/ a s ) ~ 3. The n g tion. 

Two experiments, HI and ZEUS, have reported 

their search for QCD-instanton induced processes in 

DIS at HERA HEl. While an excess of events with 

*Email: liuls@iopp.ccnu.edu.cn instanton-like topology has been observed, it cannot 
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be claimed significantly given the uncertainty of sim- 
ulation and theory. Upper limits on the cross-section 
for instanton induced processes are set dependent on 
the kinematic domain considered. 

One of the chief problems in instanton-searching is 
to identify the IFS, current jet and proton remnant. 
Only when the IFS and jet are identified precisely 
the reconstruction of the kinematic variables such as 
qi2 _ _ gr /2 j x i _ Q'2/(2g ■ q') can be good and the 
discriminating variables for instanton-searching can 
be reliable. 

At partonic level the separation of different parts 
is clear, cf. Fig. 1. The mother parton of current jet 
is the quark q" . The "mother" of IFS is the quark 
q' and gluon g in the presence of I. The best iden- 
tification of IFS and jet would be to minimize the 
differences between the energies and masses of IFS 
and current jet with those of their "mother" - 
q' + g and q" . However, this could be achieved only 
in Monte Carlo study but not in real experiment, 
because the 4-momenta at partonic level are immea- 
surable. 

Monte Carlo simulation has been applied to the 
problem in consideration by some authors 0, ITol 
111 IH El EH Their main goal was concen- 
trated on finding some methods that can be used in 
experiments, and the information at partonic level 
has been used only for evaluating the methods. 

For jet finding, both the cone algorithm 
EI ^2 an d the fc t -cluster algorithm |13L IbiL Il5j | have 
been used. The cone size has been chosen to op- 
timize the reconstruction of Q' 2 [T^j. The IFS is 
identified by a pseudo-rapidity band with width 1.1, 
where the width is chosen for the expectation of an 
isotropic IFS in its rest frame. While the k t algo- 
rithm has not done any optimization and the identi- 
fication of IFS was through quadrant [IJJ according 
to its proper position. Although these methods play 
an important role in the present searching strategy, 
the optimization needs improving. 

The aim of the present paper is to re-consider the 
Monte Carlo study of the problem and identify IFS 
through directly optimizing its energy and mass as 
well as the energy of current jet. The resulting IFS 
is of high precision and can be used in the theoretical 
study of instanton physics. The method used here 
is not directly applicable to real experiments but 
can provide instructions for the IFS identification in 
experimental study. 

The layout of the paper is as follows. In Sec. II 
some information at the partonic level of the in- 
stanton induced events will be presented. Several 
new methods for the reconstruction of IFS in Monte 
Carlo simulation will be proposed and compared in 
Sec. III. A good method is found which can recon- 
struct well the energies of current jet and instanton 
as well as the mass of the latter. The isotropy prop- 
erty of the obtained IFS will be discussed in Sec. IV 
in comparison with that of the current jet. Sec. V 
is a summary. 



II. BASIC DISTRIBUTIONS WITH 
PARTONIC INFORMATION 

Our study is based on the Monte Carlo genera- 
tor QCDINS [H E3 for instanton induced events. 
It is a MC package to simulate quark-gluon fusion 
process q' + g X in the background of instan- 
ton. It acts as a hard subprocess generator embed- 
ded in the HERWIG HI HU program. The sub- 
sequent fragmentation and hadronization are han- 
dled by HERWIG. In the following the default pa- 
rameters of the QCDINS 2.0 version are used, i.e. 
x' > 0.35, Q' 2 > 113 GeV 2 and the number of fla- 
vors is set to be n/ = 3. 

Final states in instanton induced DIS consist of 
4 subsystems: 1- scattered lepton, 2- hadrons from 
current quark, 3- hadrons created from the fusion 
of quark and gluon in the background of instanton 
{instanton final state for short), 4- proton remnant 
group of particles, cf. Fig. 1. The reconstruction is 
based on the assumption that the color forces among 
the current jet, the instanton part and the proton 
remnant are weak so that each part can be approx- 
imately separated and the momentum of hadronic 
final state of each part is approximately equal to 
those of the corresponding part at partonic level. 

In HERWIG the scattered electron is easily iden- 
tified by particle ID and status code. We discard it 
first. The particles left are the mixture of current 
jet (C for short), instanton final state (I for short) 
and proton remnant (R for short). 

After rejecting the scattered lepton, all the ob- 
jects in the hadronic final state are boosted to the 
hadronic center-of-mass frame (hem). This frame 
is defined by q + P = 0, where q and P are the 
3-momenta of the exchanged photon and proton, re- 
spectively. The positive Z-direction of hem frame is 
defined in the direction of photon momentum. In the 
hem frame, in order to see the jet signal, we remove 
the 4> angle of current quark for every final state par- 
ticles in every event. That means letting the (j> angle 
of current quark equal and hence the jet fragments 
are around = 0. The 2D distributions of 9 and (j> 
for all final state particles are plotted in Fig. 2. The 
ID distribution of 9 for all the final state particles is 
plotted in Fig. 3 (a). Also plotted in Fig.3 (b)-(d) is 
the partonic information of C, I and R. The angle 9 
is measured against the direction of photon. 

III. THREE METHODS FOR IFS 
IDENTIFICATION 

Method 1: p z -sorting 

Fig. 3(a) illustrates that the final state parti- 
cles distribute mainly around two back-to-back di- 
rections — —9 = and tt. Around 9 = are C and 
the main part of I, around 9 = ir are R and a little I. 
The method to discard R in ZEUS experiment 
is to do a 9 = ^ cut, taking particles with Q hcm > 1 
as proton remnant. In fact, as can be seen from Fig. 
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FIG. 2: (Color online) 2D plot in hem after removing 
the <f> angle of current quark. (a)(b) — cA vs. 9 in con- 
tour and surface plots, respectively; (c)(d) — cj> vs. r\ (r\ 
is defined as — log(tan|)) in contour and surface plots, 
respectively. 
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FIG. 3: The distribution of 9 defined in hem frame, 
(a) for all final state particles; (b) (c) (d) are calculated 
from the momenta of current quark, instanton part and 
proton remnant, respectively, at partonic level. 
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FIG. 4: The energy reconstruction errors of proton 
remnant by doing two cuts in 9 in hcm(a) 9 cut = 5, (b) 
The legends show the Mean and RMS. 
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FIG. 5: 9 distribution in cm2, (a) for final state particles 
(C+I), (b) for instanton before hadronization. 
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FIG. 6: The reconstruction errors for jet energy A£j ot , 
instanton energy AEi and instanton mass AMi in p z - 
sorting method. Bottom panels are the results of a cut 
on AE < 10% of top panels. 



3(a), the peak around 8 = tt is narrower than that 
around 8 = 0. So, t9 cut should be in a place on the 
right of j. We find that (9 cut = ^ is better than 
6 cui = || , which results in a better reconstruction of 
the energy of proton remnant, cf. Fig. 4. 

We define the reconstruction error of a variable 
Y as the difference between the reconstructed value 
and the true value before hadronization, i.e. 



AY 



Y„ 



Yn 



Yn 



x 100%, 



(2) 



where Yr OC represents the reconstructed value, lo 
represents the true value before hadronization. In 
contrast to the case of 8 cut = -| (see Fig. 4(a)), the 
energy reconstruction error for R after f9 cut = (see 
Fig. 4(b)) is focused at and has a smaller RMS. In 
other words, the 9 cu t = ^f- method reconstructs the 
energy of R more precisely. 

After discarding R by doing 6> cu t = the left 
particles are the mixture of C and I. Boost (C+I) to 
their c.m.s. frame (cm2 for short). Take the momen- 
tum direction of current quark, i. e. the direction of 
the jet axis, as the Z axis of cm2. The 8 distribu- 
tion of final state particles (C+I) in cm2 is shown in 
Fig. 5(a) and that of instanton before hadronization 
in Fig. 5(b). The two peaks in Fig. 5(a) represent 
current jet and IFS respectively. Since the Z axis is 
in the direction of jet axis, the peak around 8 = 
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FIG. 7: The 6 distribution in cm2 of identified C and 
I by p 2 -sorting method. 



represents the current jet particles, while the other 
peak, around 8 = 7T, represents the IFS. 

Renumber these n particles by their p z , sorting 
them according to p z i > p Z 2 > ■ ■ ■ > p zn . Ac- 
cumulating the energy from particle 1 to particle 
k gives Ek = z\ + £2 + ■■• + £fc (£j represents the 
energy of the i-th particle); simultaneously accu- 
mulating the energy from particle k + 1 to n gives 
E' k = Sk+i + £fc+2 + ••• + £«• Taking particles from 1 
to k as current jet, from k + 1 to n as IFS, the en- 
ergy reconstruction errors for jet and instanton are 

AE jct = x 100% and AE l = x 100%, 

respectively, where Eq is the energy of current par- 
ton, Ei is the energy of the quark and gluon included 
in IFS at partonic level. The value of the parameter 
k is chosen to make AE = 0.4 x |A£j ot | + 0.6 x \AEi\ 
minimum. 

The reconstruction errors for the energy of jet, the 
energy of instanton and the mass of instanton are 
shown in Fig's. 6 (a), (b) and (c), respectively. It 
can be seen that large error events occur with some 
probability and the mass of instanton has not been 
reconstructed well. The criterion on which we judge 
for a good reconstruction quality is based on that the 
reconstruction error centralizes around and has a 
narrow RMS. After doing a cut AE < 10%, both are 
improved — — large error events have been cut out 
and the mass of instanton is reconstructed better at 
the expense of throwing away 67% of the events, cf. 
Fig's. 6 (d), (e), (f). 

Up to now, the current jet and IFS are separeated. 
The 8 distribution of each part in cm2 are shown in 
Fig. 7. Jet particles centralize around 8 = 0, IFS 
centralizes around 8 = tt. At intermediate 8 the two 
parts overlap. They are separated by their p z . 

Method 2: 2D-cut 

Based on the 2D plot shown in Fig. 2 and en- 
lightened by the quadrant method (or quadrate-cut 
method) used in [14|, an improved polygonal cut is 
used in ?7-0-plane to isolate each part. Consider- 
ing that the relative position of C and I in 77-^-plane 
strongly depends on the virtuality Q 2 of photon, the 
polygonal cut is running with Q 2 , as shown in Fig. 8. 
Using these cuts each part is identified. The recon- 
struction errors are shown in the upper panels of 
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FIG. 8: (Color online) Running of the polygonal cuts 
(dashed lines) with Q 2 . The intervals of Q 2 for (a), 
(b), (c), (d), (e), (f) are [113,200], [200,400], [400,600], 
[600,800], [800,1000], [1000, 00] GeV 2 , respectively. 
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FIG. 9: (Color online) (a) (b) are 2D plots r vs. 6 in 
hem, surface and contour plots, respectively; (c) (d) are 
ID distribution for r and 6 respectively. All the variables 
are calculated in hem. 

Fig. 10. 

Method 3: r-sorting in 2D plane 

Comparing the two methods mentioned above it 
can be seen that the p z -sorting method gives better 
reconstruction. However, it disregards totally the 
information contained in the 0-0-plane, cf. Fig. 2, 
which has been taken into account in the 2D-cut 
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FIG. 10: The reconstruction errors for jet energy AE- sct , 
instanton energy AEi and instanton mass AMi in 2D- 
cut method (upper panels). Shown in the bottom panels 
are the results of r-sorting method. 



method. Therefore, we try to combine the two meth- 
ods. 

Let us define a distance in 0-0-planc, 




(3) 



where (0o,(j) = 0) is the position of current quark. 
This variable measures how far every final state par- 
ticle is from the jet axis. Choosing an appropriate 
value for tq the particles with r(9, <j>) < ro are at- 
tributed to jet. The distribution of r(0, <fr) is plotted 
in Fig. 9. 

First discard R by ID cut 9 hcm > Renumber 
the left n particles by their r, let ri < r<i < ■ ■ ■ < r n . 
Accumulating the energy from particle 1 to par- 
ticle k gives Ek = £\ + £2 + ••• + £k, simultane- 
ously accumulating energy from k + 1 to n gives 
E' k = £ k+ i + £ k+2 + ... + £„. Taking particles from 
1 to k as current jet, from k + 1 to n as IFS. The 
value of r is chosen to optimize the energy recon- 
struction, i.e. make AE = 0.4x |A£j et |+0.6x \AEi\ 
the minimum. The reconstruction errors are shown 
in Fig's. 10 (d), (e) and (f). 

It is obvious that r-sorting in 2D plane gives the 
best reconstruction among the 3 methods described 
above. 



IV. A COMPARISON OF THE ISOTROPY 
DEGREE OF IFS AND JET 

IFS are expected to be isotropically distributed 
in its rest frame. Jet fragments are anisotropy, be- 
ing axially symmetric. To quantify the degree of 
isotropy we use the sphericity. 

Sphericity is a measure of how isotropically a col- 
lection of 4-momenta is distributed in three dimen- 
sions. A normalized momentum tensor is calculated 
from the momenta of the particles in consideration, 



ga.fi = 



Epa p£ 
i i i 



(4) 



where a, /3=1, 2, 3 correspond to the x, y and z 
components. By standard diagonalization of 5 a,/3 
one find three eigenvalues Ai > A 2 > A3, with Ai + 
A2 + A3 = 1. The sphericity of the event is then 
defined as 
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(A 2 + A 3 ). 



(5) 



An ideally isotropic system with infinite multiplicity 
has S = 1, while a perfectly anisotropic system has 
5 = 0. 

For a comparison, the sphericity of IFS and cur- 
rent jet in their own rest frames identified by the 
above-mentioned r-sorting method are calculated 
and the corresponding sphericity distributions are 
shown in Fig's. 11 (a) and (b), respectively. The 
mean values are 

5i FS = 0.53 ±0.18, Sjet = 0.20 ±0.17. (6) 

Since the value of sphericity depends on multi- 
plicity, we constructed an ideal sample, which has 
the same multiplicity as the IFS sample but has 
ideally isotropic momentum distribution for every 
event, i.e. using the thermalized momentum distri- 
bution, Eq. (7) below, for each component. The re- 
sulting sphericity is drawn in Fig's. 11 (a) as dashed 
line. The corresponding Mean is S idca = 0.78 ±0.09. 
The mean sphericity for IFS nearly equal to the ide- 
ally isotropic value means that IFS is approximately 
isotropic in its rest frame. The jet is destined to get 
low value of sphericity for its axial symmetry. 

In view of the approximate isotropy of IFS we 
plot the distribution of the three momentum compo- 
nents p x ,py,p z of IFS in its cm. frame as shown in 
Fig. 12 (a). They coincide approximately. The dis- 
tribution of the average pi of these three components 
is shown in Fig. 12 (b). 

For a thermalized system the momentum compo- 
nent distribution is 



1 



^2nmk B T 



e -p>/Vrnk B T) dpu 



i = x,y,z. 
(7) 

Fitting the pi distribution shown in Fig. 12 (b) to 
this formula we get the "temperature" of IFS as 

T = 863.86 MeV for m = , 
= 244.22 MeV for m = m K , 
= 128.50 MeV for m = mP . (8) 

This "temperature" is actually not a thermal one, 
since at such high temperature there will be decon- 
finement and no hadron can survive. 

The "temperature" T in Eq's. (8) is a parameter 
characterizing the phase-space distribution of par- 
ticles in IFS. Its high value means that the quark- 
gluon fusion in the presence of instanton is a very 
hard process. 
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FIG. 11: Sphericity for IFS (a) and current jet (b) in 
their respective rest frame. Dashed line in (a) shows the 
sphericity for a ideally isotropic model with the same 
multiplicity as the IFS sample. 
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FIG. 12: (a) Momentum-component distribution of IFS 
in its rest frame; (b) the distribution of the mean value of 
the three components, fitted to thermalized momentum 
distribution. 



The identification of instanton final state, current 
jet and proton remnant in instanton induced deep 
inelastic scattering is studied using Monte Carlo sim- 
ulation. Different methods - - p z -sorting, 2D-cut 
and r-sorting, for the reconstruction of the quark- 
uon fusion product in the background of instanton 
(IFS) and current jet are tried and compared using 
QCDINS Monte Carlo event generator. A method to 
optimize the energy reconstruction is applied, which 
can reconstruct well the energies of current jet and 
instanton as well as the mass of the latter. It is found 
that r-sorting method gives the best reconstruction. 

The sphericities of IFS and current jet identified 
by the r-sorting method are calculated. The high 
\alue of sphericity for IFS means that IFS is ap- 
Pl i roximately isotropic in its rest frame. 

The momentum distribution of IFS approximately 
riimics a thermalized distribution. The "tempera- 
ture" fitted is rather high, which can be taken as 
a charateristic parameter for measuring the "hard- 
r ess" of the quark gluon fusion process in the pres- 
ence of instanton. 
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